Density of States of Disordered Two-Dimensional Crystals with Half-Filled Band 
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A diagrammatic method is applied to study the effects of commensurability in two-dimensional 
disordered crystalline metals by using the particle-hole symmetry with respect to the nesting vector 
Po = {±-, ^} for a half-filled electronic band. The density of electronic states (DoS) is shown to 
have nontrivial quantum corrections due to both nesting and elastic impurity scattering processes, 
as a result the van Hove singularity is preserved in the center of the band. However, the energy 
dependence of the DoS is strongly changed. A small offset from the middle of the band gives rise to 
disappearence of quantum corrections to the DoS . 
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i" 1 ' Great advance has been made in the theory of dis- 
ordered metals after the pioneering work of Abrahams 
et al. |l| . According to this paper all electronic states in 
one- and two-dimensional (Id and 2d) disordered systems 
are_ localized irrespective of the degree of randomness, 
Electron-electron correlations in disordered met- 
lave been shown to result in nontrivial corrections to 
the density of electronic states, (DoS), and conductivity, 
The corrections to the conductivity are similar to 
Q ' the localization corrections obtained for a noninteracting 
O | electron gas, and quantum corrections to the DoS were 
shown to reduce it near the Fermi level. 
' Disordered metals in all abovementioned papers are 
^ . modeled as a free electron gas moving in the random 
field of rigid impurities. However, at low concentrations 
ON , of impurities the crystal usually has a periodical struc- 
ture and the impurity atoms in most cases substitute the 
host atoms of the lattice. Then, the effects of commen- 
surability of the electron wavelength, A , and the lattice 
| constant, a, become essential in the scattering processes. 
. The commensurability is known to exist at all 'rational' 
points of the electron band, however it appears to be im- 
£h , portant for a half-filled band. 

I 1 In this Letter we present our study concerning the ef- 
fects of weak disorder on the electronic DoS of a 2d crys- 
talline metal. The Hamiltonian of the model can be writ- 
ten as H = H a + V(r), where H a is the Hamiltonian of 
noninteracting electrons in the perfect square lattice with 
nearest- neighbor hopping and V(r) = ^ . U(r — Rj) is 
the impurity potential with Ri being the positional vec- 
!_h ' tor of an impurity randomly located on the i — th lattice 
site. 

The one-particle DoS of the regular lattice can be ex- 
pressed as pf ] = {2 ^ R)d J |v ^ fc)| , where dS is the ele- 
ment of an isoenergetical surface in d-dimensional space. 
It can be shown that p^f\e) has a van Hove singularity 
at the points where the group velocity of the electron 
wave packet Vk = Ve(fe) vanishes, H. For a three- 

(3) 

dimensional (3d) regular lattice Pq (e) has integrable sin- 
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gularities. 

For a pure 2d lattice with nearest-neighbor hopping 
the van Hove singularity has logarithmic character. The 
simplest electron spectrum for a 2d square lattice can be 
written in the tight-binding approximation as 

e(fe) = t[2 — cos(k x a) — cos(k y a)]; (1) 

wheve,k XiV = ^^-n XtV with < n xy < and 

only electron tunneling between nearest-neighboring sites 
with the tunneling integral t is involved. The bandwidth 
is W = At; and for a half-filled band case the Fermi 
energy becomes cf =2t. The DoS of a 2d square lattice 
with nearest-neighbor hopping is expressed by the elliptic 
integral of the first kind, an asymptotic expression of 
which has a logarithmic singularity in the middle of the 
energy band as: 

r4,t 2 -\e\ 



where e is an electron energy measured from the Fermi 
level e = e — 2t. (Hereafter the tilde on e will be dropped). 

The DoS of a noninteracting electron gas moving in the 
random field of impurities has no essential singularities 
near the Fermi surface. Inclusion of even short-range 
correlations in the 2d disordered metal give s rise to a 
decreasing DoS near the Fermi level, |p|,|fi|Jlo|. The DoS 
of a 2d disordered crystal with substitutional impurities 
turns out to have a singularity near the middle of the 
band even for the noninteracting electron gas. 

Notice that the commensurability effect for a Id dis- 
ordered crystal near the middle of the band has been 
studied by many authors, 11-^]. Dyson first pointed 



out JLIJ that the phonons' DoS of a Id disordered chain 
has a singularity as p' 1 ' (e) oc — |e| 1 ln~ 3 |e| near the mid- 
dle of the band. Later an analogous singularity has been 
found in the electronic DoS of many Id models, [[j~2] H . 
However, there exist a few computational studies of the 
DoS of a 2d disordered crystal, |20|-|2^| . By studying the 
averaged Green's function for a disordered system with 
n orbitals per site, the expansion coefficients in power of 
1/n for d > 2 were shown in [ p4| to diverge for energies 
approaching the band center. This fact was interpreted 
in [^4| as a existence of a van Hove singularity in the 
DoS. 

The technical difficulties in analytical calculations of 
physical parameters are connected with the cosine en- 
ergy spectrum and, furthermore, with an absence of per- 
turbativc parameters for half-filling. It is appropriate to 
notice that the recent attempt to explain both the lin- 
ear resistivity and the nearly temperature independence 
of the thermopower in the cuprate superconductors |25| l 
is based on the existence of a van Hove singularity in 
the DoS of these materials. To study the problem an- 
alytically, the idealized model [^6|,^7j e(fc) = k x k v for 
the energy spectrum, which gives a logarithmic DoS at 
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half-filling, was used. There exists also the nested Fermi 
Liquid scattering approach |i8|]2g| l to study the suscepti- 
bility of high-T c superconductors. 

As it is known, the Fermi surface of an infinite 2d lat- 
tice with nearest-neighbor hopping is changed with band- 
filling and it is flat for the half-filling. In this case there 
exists a nesting vector P = {±f, f } that maps an en- 
tire section of the Fermi surface onto another, i.e. the 
Fermi surface is perfectly nested for the half-filled band 
case. There, the following particle-hole symmetry of the 
electron dispersion with respect to the vector P Q for a 
half- filled band holds: 

e(p + P )-e F = -[e(p)-e F ] (3) 

The one-electron DoS can be calculated according to 
the following expression: 

p(e) = ~Tm J -^G R {p,e) (4) 

where Gn(p,e) is the retarded Green's function. 

The new class of diagrams which gives an essential con- 
tribution to the DoS is drawn in Fig.|l|. Thin solid and 
dashed lines in Fig.l correspond to the 'bare' retarded 
Green's functions Gr°(p, e) and Gr°(p + P Q , e) . 




CV(k,e) 
^ N -p + k + P c 



Ar(k,e) 




Ar(k,e) 




b) (b') " c) (c 5 ) 

FIG. 1. Self-energy parts which give a contribution to the 
Green's function. Impurity scattering is denoted by a dotted 
line with a cross. C w , CV and D^, are cooperon and 
diffusion blocks, respectively. The diagrams b') and c') are 
the symmetrical conjugate to those shown by 6) and c) with 
respect to the additional impurity line. 



Averaging over impurity realization is performed ac- 
cording to the crossed diagram technique in the Born 
approximation described in f3(i|| , i.e. the white-noise im- 
purity potential is used which corresponds to the off- 
diagonal disorder in the problem. The main diagrams 
are selected according to the condition of e F r 3> 1 or 
kpl 3> 1 with I = vpT and r being an elastic scattering 
time, which means I 3> A. 

The scattering processes in the problem are described 
by two relaxation times r G and tv, (Fig. 2). t cor- 
responds to the normal scattering process contribution 
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to which gives the diagram shown in Fig. 2a and — = 

8^ / m\ u ( s ^> & Here ' c ™p and u arc thc ° im - 

purity concentration and potential, respectively. The di- 
agrams which give a contribution to Green's functions 
contain new impurity vertices, shown in Figs. 2b, c. These 
vertices characterize simultaneous reflections of electrons 
on the Brillouin zone boundary in the process of scatter- 
ing on an impurity. The following expression represents 
these vertices: ^ = / J^\U(P , S)\ 2 . It can be 
seen from Figs. 2b, c that the total momentum, generally 
speaking, is not conserved for these impurity vertices. 
Therefore they represent an U mklapp process. 

jj - p + k_ - P' + M- P p±k pi j- k + p„ 



P> c P,C P + Po P' P + Po p' 

a) b) c) 

FIG. 2. 

Bare retarded and advanced Green's functions 
G° R A (p, e) = [e— (e(p) — ep) + ^signe] -1 contain a total 
relaxation time, which is equal to — = — + — . 

The maximally crossed diagrams are redrawn in Fig. la 
so that the self-energy part contains the particle-particle 
propagator CV(fc, e)which we call 7r-cooperon, (Fig. 3a). 
In contrast to the cooperon for an isotropic system 
C„(k, e) has a diffusion pole at a small total energy and 
a large total momentum ~ P Q . In the limit of small e 
and k , satisfying the conditions |e|r <C 1 and kl -C 1, 
C 7r (k, e) has the following form 

r (k a — c imp \u(p )\ 2 

W>*)- i_4ir|e| + (fc/)2-(r/r w )2' (&j 

which has a diffusion pole for r w — > r. Notice that e is 
an energy measured from the Fermi level. 

The self-energy parts in Figs.l b(b') and c(c') contain a 
particle- hole propagator D 7r (k,e) which can be written 
for |e|r <C 1 and kl <C 1 as 

n (k A — C imp \U{P )\ 2 

7r- diffuson Djj(k, e) has a diffusion pole in the particle- 
hole channel at a small total energy and large momenta 
difference ~ P Q for t —> r n . 

According to the Dyson equation, the retarded Green's 
function Gr(j>, e) is expressed in the following form: 

1 °° 



n=0 



(J) 



where, E(p, e) = Yli— a ^i(P> e )' and eacn 01 Si(p, e) cor- 
responds to the one of self-energy parts a — c' in Fig.l, 
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respectively. By summing all diagrams in Fig.l the self- 
energy part S(p, e) is reduced to the form: 



S(P ' £)= / ^ {Cw(fe ' c)_ ^ (1_ ^ )Ar(fe ' 
G° fl (-p + k + P ,e); 



e)} 



(8) 



-p + k + P ,e P>,e 

C ff (k, e) = | 



p, e -pi + k + i*T 



a(k,e) 



p + k + P P', e 

-> 



* + 



p,e p/-k + P ,e 



§S3 



FIG. 3. Ladder diagrams for the 7T- cooperons C\(k,e) 
and CV(k,e), and the ir- diffusons _D^(k,e) and Z),r(k, e). 
C n (k,e) and -D^(k, e) are expressed by C n (k, e) and D„-(k, e) 
as C T (k,e) = -^C w (k,e) and D T (k,e) = -^fl,(k,e) for 
small k and e. 

Substituting Eqs.(7) and (8) into Eq.(4) the following 
expression for p(e) can be obtained after integrating the 
bare Green's functions over p and taking a sum over n : 



p(e)=p^{l-Re 



4T 2 a(e) 



V / l + 4T 2 a(e)(l + y/1 + 4r 2 a(e)) 



} 

(9) 



where 



a(e) 



/ 7^{^^)-?(i-r)A ; !fc.i}: (10) 



(2tt) 2 



The obtained expression for the DoS contains high 
order logarithmically divergent contributions. So, infi- 
nite order impurity blocks are summed up according to 
Eqs.(7)-(8). 

The nesting processes are weakened away from the 
middle of the band. In this case the normal scattering 
processes become more probable and t q <C . As a 
result the diffusion poles in the propagators C n (k, e) and 
D^{k,e) disappear. Therefore a small offset from the 
middle of the band gives rise to the disappearence of the 
quantum corrections to the DoS and p(e) = p . When 
filling reaches the center of the band the Fermi surface 
becomes flat. In this case the scattering with nesting 
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seems to have preference, i.e. tv <C t q and r = t v . The 
impurity blocks C 7r (k.e) and D^ik^e) have a diffusion 
pole under this condition and we get from Eq.(lO): 

«( e ) = o 1 s ln(--p— j-r) (11) 
87T6 ft-k 4z7v|e| 

Substitution this expression for a(e) into Eq.(9) gives 
the following expression for the DoS: 

p( e ) = pi 2) W{i--==^ 7 =—}; (12) 

y/l + a (e)[l + y/l + a (e)] 

where,a Q (e) = 2 ^ FT ^ \n(j^). Near the vicinity of 
the Fermi level Eq.(12) is approximated as 

^ = ^W^^4^R)]" 1/a (13) 

p£\e) in Eqs.(12) and (13) is the DoS of a pure 2d 
square crystal, which has the van Hove singularity ex- 
pressed by Eq.(2). Elastic scatterings in the crystal with 
substitutional impurities preserve the central peak of the 
DoS, however the energy dependence of p(e) is changed 
from logarithmic dependence to the square root of the 
logarithm in the close vicinity of the middle of the band: 

p(e)^-^(2ne F T^ 1/2 \n^ 2 (-\-) as |e| - 
(na) €f 4T 7rl e l 

(14) 

The peak also becomes narrower than that of the van 
Hove one due to impurity scattering. 

It is worth to compare here the result presented in 
the Letter with that obtained for a Id lattice with 
half filled energy band, containing off-diagonal disorder, 
@|2)-[l3Jl§-g. The DoS of the pure Id lattice is a 
smooth function of the energy within the band and it has 
a singularity only at the boundary of the energy band. 
Therefore the Dyson peak of the DoS in the middle of the 
half-filled band of the Id disordered lattice is a result of 
strong Bragg reflection of the electrons in the process of 
scattering on impurities with a consequent interference of 
scattered waves. The calculation of density-density and 
current-current correlators for Id disordered chains shows 
that the localization length diverges and static conduc- 
tivity saturates to a constant value at T = with the 
approaching of half filling. In contrast to the Id case the 
singular enhancement of the DoS in the 2d disordered lat- 
tice is due to a van Hove singularity in the bare DoS. Ef- 
fect of impurity tends to decrease the DoS. The obtained 
result is also confirmed by the computational study of 
the DoS of a 2d square lattice with off-diagonal disorder 
[ 23 1 . Numerical analysis of the participation numbers in 
[23 1 shows that localization becoms stronger close to the 
band center. 

Notice that a nested Fermi surface with a nesting vec- 
tor P* = £P„ (£ < 1) arise for a tight-binding mod- 
els which include both nearest neighbor and next near- 
est neighbor hopping terms and correspond to a non-half 
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filled band. In this case the van Hove singularity in the 
DoS of pure system disappears and the DoS vanishes on 
the Fermi surface. 

In conclusion, we proposed a new diagrammatic ap- 
proach to study nesting effects on the DoS of 2d disor- 
dered square lattice. Calculations were performed in the 
diffusion approximation. However, unlike the "conven- 
tional" localization theory, the infinite order logarithmi- 
cally divergent diagrams were summed up here according 
to the Dyson equation. The quantum corrections tend to 
decrease the DoS in the whole energy interval. As a re- 
sult, the van Hove singularity in the middle of the band 
is preserved, however, its energy dependence is changed 
from logarithmic to square root of logarithm in the close 
vicinity of the band center. The calculations are valid 
for a half-filled energy band. A small shift from the band 
center results in the vanishing of the quantum corrections 
to the DoS. 

The authors are indebted to a referee who pointed 
out the difference between impurity effects in Id and 2d 
systems. One of the authors (E.P.N.) thanks A.Erzan 
and all members of the Physical Department of Istan- 
bul Technical University, where the part of this work was 
performed, for long time hospitality. E.P.N, also thanks 
H.Feldmann for discussion. 
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